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Abstract — This paper maximizes the achievable throughput 
of a relay-assisted wirelessly powered communications system, 
where an energy constrained source, helped hy an energy 
constrained relay and hoth powered hy a dedicated power heacon 
(PB), communicates with a destination. Considering the time 
splitting approach, the source and relay first harvest energy 
from the PB, which is equipped with multiple antennas, and 
then transmits information to the destination. Simple closed-form 
expressions are derived for the optimal PB energy heamforming 
vector and time split for energy harvesting and information 
transmission. Numerical results and simulations demonstrate the 
superior performance compared with some intuitive benchmark 
beamforming scheme. Also, it is found that placing the relay at 
the middle of the source-destination path is no longer optimal. 


I. Introduction 

Gigabit wireless access will be a reality in fifth generation 
(5G) wireless systems, with a series of breakthroughs, such 
as massive multiple-input and multiple-output (MIMO), full- 
duplex communication and small cell architectures, which 
in turn, has fueled a number of emerging wireless services 
such as mobile gaming, mobile TV and mobile Internet. 
With the proliferation of smartphones and tablets, one of 
the most critical issues affecting the user experience is the 
limited operation lifetime of mobile devices due to hnite 
battery capacity. Motivated by this, radio-frequency (RF) based 
energy harvesting technique has received substantial research 
interests in recent years ID, El- Empowered with the RF 
energy harvesting capability, it is possible to virtually provide 
perpetual energy supply to mobile devices, eliminating the 
need to plug into the power grid for battery recharging. 

RF energy harvesting, combined with information transfer, 
has resulted in a new emerging topic, namely, simultaneous 
wireless information and power transfer (SWIFT), which has 
attracted signihcant attention from academia and recently 
from industry. Thus far, various aspects of SWIFT systems 
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have been investigated, including information theoretical limits 
a, a, practical architectures a-Q, effect of imperfect 
channel state information a, OFDM-based SWIFT systems 
a, IfTOl and relay-assisted SWIFT systems IfTTI - lfTTl . It is 
worth pointing out here that all these prior works assume 
that the mobile devices harvest energy from ambient RF 
signals. However, even if this may be viable for low power 
devices such as sensors, it is in general infeasible to power 
larger devices such as smartphones, tablets and laptops na. 
Responding to this key limitation, the authors in ifT^ proposed 
a novel network architecture, where cellular base stations are 
underlaid by dedicated power beacons (FB), which can be 
used to supply energy to mobile devices through microwave 
power transfer. Since the FBs do not require any backhaul link, 
the associated cost of FB deployment is much lower, hence, 
dense deployment of FBs to ensure network coverage for a 
wide range of mobile devices is feasible. 

In this paper, we consider a dual-hop decode-and-forward 
(DF) relaying system, where both the source and relay are 
powered by the dedicated FB. To improve the energy transfer 
efficiency, the multiple antenna enabled FB performs energy 
beamforming. It is assumed that wireless power transfer is 
performed over the same frequency as information transfer, 
and the time-switching protocol 13 is adopted. As such, 
the single antenna at the source and relay switches between 
two hardware chains, one for energy harvesting and one 
for information transmission/reception. Therefore, the entire 
communication block consists of two different phases, i.e., 
energy harvesting phase, where the source and relay harvest 
the energy from the FB, and information transfer phase, where 
the relay assists the information transmission between the 
source and destination. 

The main contribution of this paper is the derivation of 
simple closed-form expressions for the optimal time split 
between the energy harvesting and information transfer phase, 
as well as the optimal energy beamforming vectors at the 
FB, which maximize the system’s throughput. Simulation 
results demonstrate that the optimal solution substantially 
outperforms the intuitive benchmark beamformer. In addition, 
it is revealed that placing the relay in the middle of the 
source and destination path is no longer optimal, instead, the 
distances to the FB should also be taken into consideration 
when optimizing the relay position. 

Notations: (•)* denotes the complex conjugate, de¬ 
notes the matrix transpose, and ( denotes the Hermitian 
transpose. I is the identity matrix of appropriate size. Hx = 
X(X^X)“^X^ is the orthogonal projection onto the column 
space of X, and H^ = I — Hx is the orthogonal projection 
onto the orthogonal complement of the column space of X. 
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II. System Model 

We consider a communication system where the source S 
communicates with destination D with the assistance of the 
relay R, as depicted in Fig.|3 We assume that both S and R are 
energy constrained, hence rely on the external energy charging 
via wireless power transfer from a dedicated PB. All three 
communication nodes are equipped with a single antenna, 
while the PB is equipped with N antennas. Full channel state 
information (CSI) of the PB and source and relay links is 
assumed at the PB. In practice, the channel can be estimated by 
overhearing the pilot send by the source and relay. In addition, 
the channel magnitudes of the two hop information links are 
assumed to be known at the PB, which, for instance, could be 
obtained through feedback from the relay. 



In the first half of the second phase, S transmits information 
to R using the energy harvested in the first phase. Hence, the 
received signal at R is given by 


2E.. 


Vr = 


(1 - r)Tdf 


/iSO ~\~ Tlr, 


(3) 


where denotes the distance between S and R, fi is the 
channel coefficient, sq is the information symbol with unit 
energy. 

Upon receiving the source signal, the relay first decodes the 
source symbol and then forwards it to the destination using 
the energy harvested during the first phaseQ as such the signal 
at D can be expressed as 


Vd 


I 

V (1 - 


+ nd, 


(4) 


where ^4 denotes the distance between R and D, /2 is the 
channel coefficient, Sr is the relay signal with unit energy, 
and Ud is a sample of the AWGN with variance Nq. 

Therefore, the effective end-to-end SNR at the destination 
can be written as 


2Tr]P . [ |w^hip|/i|2 |w^h2p|/2p 


(5) 


Fig. 1; System model 


III. Throughput Optimization 


The entire communication consists of two different phases, 
namely, energy harvesting and information transmission 
phase. Assuming a block time of T, during the first phase 
of duration tT, where 0 < r < 1, S and R harvest energy 
from the PB. The remaining time of duration (1 — t)T is 
equally partitioned into two parts, during the first half period, 
S transmits information to R and during the second half, R 
forwards the information to D. 

During the energy harvesting phase, the received signal at 
S and R can be expressed, respectively, as 

vl = \JPdf + ris and y® = -f 


The achievable system’s throughput, R, can be expressed 
as 


(1 — r) 

R = ^ 


2 
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Hence, the following optimization problem is created: 


^ 1/1 


|w^h 2 | 2 |/ 2|2 


dfd^ 
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r.w 

R 

(7) 


S.t. 

0 < r < 1 

(8) 



||wf = 1. 

(9) 


where P is the transmit power at the PB, di and d 2 denote 
the distances between PB and S, PB and R, respectively. 
Furthermore, a is the path loss exponent, hi and h 2 are the 
channel vectors of size W x 1, Xe is an x 1 vector satisfying 
E {x|xe} = 1, and and Ur are the zero-mean additive white 
Gaussian noise (AWGN) samples with variance Nq. 

Since the PB is equipped with multiple antennas, energy 
beamforming could be applied to improve the efficiency of 
energy transfer, i.e., Xe = wse, where w is the beamforming 
vector with ||w|p = 1 , while Sg is the energy symbol with 
unit power. As such, the total received energy at the S and R 
at the end of the first phase can be expressed as 

f;, = p|w^hi|2prrdr“, (1) 

and 

= p|w^h 2 pPTTd^“, ( 2 ) 

respectively, where 0 < 77 < 1 is the energy conversion 
efficiency. 


At the first glance, the above problem requires the joint 
optimization of r and w, which is in general difficult. Never¬ 
theless, a close observation reveals that the special structure of 
Vl allows for a separate optimization of t and w. Specifically, 
we present the following key result. 

Proposition 1: The original optimization problem, Vl, is 
equivalent to the following: 

P2: max log, (^1 + (10) 

S.t. 0<T<1, (11) 


where Zm is defined as 


Zm = min 


Iw^hiplAp |w^h2n/2p ] 

dfd’^ ’ I ’ 


( 12 ) 


* Please note, we have ignored the processing power required by the 
transmit/receive circuitry at the relay as in Cl, im, (20) This assumption 
is justifiable since the transmission energy is the dominant source of energy 
consumption. 
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with w being the solution of the following optimization 
problem: 


: 


max min 

W 


Iw^hiplAp 

’ d^d2 I 


s.t. ||wf = l. 

Proof: Define function g{T, z) as 

g(r,A = ^^^log2 (l + ' 


2Tr]P 

(1 - r)No 


(13) 

(14) 


(15) 


where 0 < r < 1 and z is a positive real number. It is easy 
to prove that function g{T, z) is an increasing function with 
respect to z. Now consider two positive real numbers zi and 
Z 2 such that Z 2 > zi, and let ti being the value of r which 
maximizes g{T, zf), i.e., p(Ti, zf) > gir, Zi) for all t, i = 1, 2. 
Then, it holds that 


9{t2,z 2) > g{Ti,z2) > g{Ti,zi), (16) 


which indicates that the maximum of g{T, z) is achieved at 
the point when z attains its maximum. Therefore, a sequential 
optimization of problems V3 and V2 yields the optimal 
solution for the original problem VI. ■ 

In the following, we investigate the optimal solutions for 
the problems V2 and VS. 

Proposition 2: The optimal r for the optimization problem 
V2 is given by 


T = 




(17) 


where W{x) is the Lambert W function IItTI . and /3 = ■ 

Proof: The proof follows from the results presented in 
ifTSl Appendix A]. ■ 

We now turn to problems VS, and we have the following 
key result: 

Theorem 1: The optimal beamforming vector w for the 
optimization problem VS can be expressed as 




W = X- 


inh.hil 


+ '^1 — X 


^2 


KM 


(18) 


where hi = —|iii=hi, ±±v — , 

s/d-id^ s/d^ ^4 

|n^ hi||, and c = and x being 


h2 = -^h2, a = Wlli^h.W, b = 


h 2 




X = 


a „ > af+P 

Va^+b^ ’ - “ ’ 

, ^-, a<c< 2 !+^ 

i/(a-c)2+f,2 ’ - a 

1 ^ ^ ^ 


(19) 


Proof: According to fT2\. Corollary 1], the optimal beam¬ 
forming vector w can be expressed as 


w = X- 




+ \/l — 


n^.h! 

KKw 


( 20 ) 


where x £ [0,1]. Now let us define gi{x) = jw^hi] = 
ax + bVT — x"^, and g2ix) = |w^h 2 | = cx. Then, the 
original optimization problem 7^1 is equivalent to maximize 


the function g{x), i.e., 

max g{x), s.t. 0 < X < 1 , ( 21 ) 

X 

where g{x) is defined as 

g{x) = min {gi{x), g2{x)) = min ^ax + b \/1 — x^, cx^ . 

It is easy to show that gi (x) is a concave function with respect 
to X, hence, its maximum can be attained by solving g[ (x) = 
0 , which gives x = -^=^= 5 , and gi{x) = \/af~+lP. 



(c) Case 3 


Fig. 2: Three different cases for the maximization of g{x). 

Now, as shown in Fig. we have three different scenarios: 

• In Case 1, we observe is that the x-label of the cross 

point of gi{x) and 52 ( 2 :), i.e., Xc, can be characterized 
by 5 i( 2 ;c) = 92 {xc), which gives x^ = If 

the slope of 52 (x) is sufficiently large, such that the cross 
point appears before 9 i{x) attains its maximum, i.e., Xc < 
X, which gives c > ^ , then, the maximum of g{x) is 

achieved at x = x = , which is also the maximum 

Va2+62 _ 

point of 5 i(x), as shown in Fig. |2(a)| 

• In Case 2, the cross point appears after 9 i{x) attains its 

2_|_r2 

maximum, i.e., a < c < 2_zt2_ then, the maximum 
of g{x) is achieved at the cross point, i.e., x = Xc = 
^^=L_, as shown in Fig. ^ 

• In Case 3, there is no cross point between pi(x) and 
52 ( 2 ;), be., 9 i{x) > 92 {x) for all x, namely, c < a, then, 
the maximum of g{x) is identical to the maximum of 
52 ( 2 :) which is achieved at the point x = 1, as shown in 

Fig.M 


IV. Numerical Results and Discussion 

In this section, numerical and simulations results are pre¬ 
sented to illustrate the impact of key system parameters on 
the system’s throughput. Without loss of generality, we set the 
energy conversion efficiency rj = 0.4, and path-loss exponent 
a = S. Please note, the throughput is obtained by averaging 
over 10 ^ independent channel realizations. 
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To demonstrate the superiority of the proposed optimal 
scheme, we compare it with an intuitive benchmark scheme 
by looking into the asymptotic large antenna regime, where 
the optimal beamforming vector becomes 

- { I/2I ht , I/ll h- ^ / /I/2P , l/iP 

ydpf llhill ||h2|| y / V df df d“d“ ■ 

( 22 ) 

The rationale behind the choice of w is that, the optimal 
beamforming vector should be a linear combination of 
and h2, hence, the key is to design the optimal weights. 
Capitalizing on the asymptotical orthogonality of and h| 
when N 00, the optimal weights can be easily obtained. 

Fig. 13 depicts the achievable throughput of the optimal 
scheme and the benchmark scheme with di = d 2 = 3m 
and d^ = di = 5m. It can be readily observed that the 
optimal scheme outperforms the benchmark scheme, and the 
performance gap is rather significant for moderate number of 
antennas N. On the other hand, when N is sufficiently large, 
i.e., N = 5000, the performance gap narrows substantially. 
This is rather expected since the benchmark scheme becomes 
asymptotically optimal. We also observe that the throughput 
improves when N increases, which is also intuitive since the 
energy transfer efficiency improves with a large size of antenna 
array. 



Fig. 3: Throughput comparison between the optimal scheme 
and benchmark scheme for different N. 

Fig. 0] examines the impact of node positions on the 
throughput performance when N = 10, di + ^2 = 20m, and 
d^ + di = 20m, as a function of di and d^, both vary from 
7 to 13. It can be observed that higher throughput is attained 
at the point di = 7 and ds = 13, a scenario where the PB is 
close to the source while the relay is close to the destination; 
or at the point di = 13 and da = 7, a scenario where the PB 
is close to the relay while the relay is close to the source. The 
above results are somehow intuitive, since the performance 
of dual-hop relaying systems is bottlenecked by the weakest 
link, hence, an optimized system shall achieve a fine balance 
between the two hops. 

Another interesting observation from Fig. 0] is that the 
symmetric setup, i.e., the point di = d2 = 10 and da = 
di = 10, does not yield the maximum throughput. This is in 
sharp contrast to the conventional dual-hop relaying systems 
where it is always desirable to put the relay node in the 
middle of the source and destination link. The reason is that, 
with the introduction of PB, in addition to the distance of 


the information transfer links, the throughput performance 
also heavily depends on the distance of the energy transfer 
links. As a matter of fact, the throughput is determined by 
d^dg as shown in the end-to-end SNR expression (|5]l. Since 
(7 X 13)“ < (10 X 10)“, it becomes obvious why the maximum 
throughput is achieved at point di = 3 and da = 13. 
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Fig. 4: Impact of node positions on the system throughput. 

Fig. |5] compares the achievable throughput of the dual¬ 
hop relaying system with that of direct transmission system 
with optimized time split and beamforming vector. As can 
be readily observed, at the SNR levels of practical interest, 
i.e., 0 dB < SNR < 50 dB, adopting the relaying structure 
improves the system throughput. Moreover, the performance 
gap is more substantial with moderate number of antennas, 
i.e., N = 10. Only if the transmit SNR is very high, direct 
transmission becomes preferred. This is rather intuitive since at 
the high SNR regime, the system is degree-of-freedom limited, 
as such the half-duplex relay operation becomes the bottle¬ 
neck as manifested through the 1/2 factor in the throughput 
expression. Since in the wirelessly powered communications 
systems, the source is likely to operate in the power-limited 
regime, adopting the relaying structure is beneficial in general. 



Fig. 5: Throughput comparison: Relaying v.s. Direct 
transmission. 


V. Conclusion 

In this paper, we have optimized the throughput of a relay- 
assisted wirelessly powered communication system. Specifi¬ 
cally, we obtained simple closed-form solutions for the PB 
energy beamforming vector as well as the optimal time split 
for the energy harvesting phase and information transmission 
phase. It was shown that the optimal solution yields signif¬ 
icant performance gain compared to the intuitive benchmark 
scheme. 
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